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ABSTRACT 

We present an analysis of ground-based optical photometry and spectroscopy, and Rossi X-ray 
Timing Explorer X-iay observations of the old novae DI Lacertae and V841 Ophiuchi. Our optical 
photometry data (obtained with the automated photometry telescope RoboScope) comprise an 
almost decade-long light curve for each star, while the contemporaneous spectroscopy and X-ray 
observations repeatedly sampled each nova during separate intervals of w 45-55 d in length. 
The long-term optical light curves of both novae reveal quasiperiodic variability on typical time 
scales of ~ 30-50 d with amplitudes of AV ~ 0.4-0.8 mag. V841 Oph also displays a long-term, 
sinusoidal modulation of its optical light on a time scale of 3.5-5 yr. The optical spectra of these 
novae display quite different characteristics from each other, with DI Lac showing narrow Balmer 
emission cores situated in broad absorption troughs while V841 Oph exhibits strong single-peaked 
Balmer, He I and He II emission lines. We find little change between spectra obtained during 
different optical brightness states. The X-ray count rates for both novae were very low (< 1.5 ct 
s~^) and there was no reliable correlation between X-ray and optical brightness. The combined 
X-ray spectrum of DI Lac is best fit by a bremsstrahlung emission model (with kT ^ 4 keV and 
A^H < 1.8 X 10^^ cm~^); the X-ray spectrum of V841 Oph is too weak to allow model fitting. We 
discuss the possible origin of variability in these old novae in terms of magnetic activity on the 
secondary star, dwarf nova type disk instabilities, and the "hibernation" scenario for cataclysmic 
variable stars. 
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1. Introduction 

Cataclysmic variables (CVs) are semi-detached 
interacting binary stars composed of a white dwarf 
(WD) primary star and a low mass (< O.5M0) 
main sequence secondary star, with typical orbital 
periods of < 1 d. The Roche-lobe-filling secondary 
star loses mass through the inner Lagrangian point 
into an accretion disk formed around the WD (for 
non- magnetic CVs). Classical novae are a sub- 
class of CV in which a thermonuclear runaway is 
triggered in a reservoir of matter that has been 
gradually accreted onto the WD. The resultant 
outburst produces a peak brightness increase of « 
6-15 mag, and releases 10'''* — lO'*^ erg (see review 



in Warner 1995, ch. 5). 

The outbursts of novae are often well-observed, 
and the long term behavior of many dwarf novae 
and novalike CVs is monitored by the American 
Association of Variable Star Observers. However, 
there have been few programs for monitoring the 
variability of novae at quiescence. One such pro- 
gram, using an automated telescope called Ro- 
boScope (see §2.1), has monitored 22 old novae 
and 42 novalikes, most of them for over 9 years 
now. RoboScope has found several kinds of un- 
usual photometric behavior. About one-third of 
the old novae have shown quasi-periodic variabil- 
ity for a year or two, interspersed with stable light 
curves. These variations do not appear stochas- 
tic as they repeat at similar periods when they 
reappear. Preliminary analysis of RoboScope light 
curves spanning several years for the old nova 
DI Lacertae suggested brightness oscillations on 
a time scale of « 35 d (Honeycutt et al. 1995), 
while the old nova V841 Ophiuchi also displayed 
prominent variability (Honeycutt et al. 1994). In 
order to determine if these oscillations could be re- 
lated to a disk instability mechanism (as operates 
in dwarf novae) or to mass transfer or magnetic 
effects on the secondary star, we attempted a de- 
tailed study of DI Lac and V841 Oph, which are 
among the most active and brightest of the sys- 
tems showing the oscillations. 

DI Lac (= Nova Lac 1910) was a moderately 
fast nova that reached a maximum brightness of 
my(max) = 4.3 mag during its outburst. V841 
Oph (= Nova Oph 1848) was a slow nova that 
reached a similar brightness, my (max) = 4.2 mag 
(Kukarkin ct al. 1971). The reddening for both 
objects has been measured from International Ul- 
traviolet Explorer (lUE) ultraviolet and ground- 
based optical spectra. The range for DI Lac is 
E{B -V)= 0.15 0.41 (Bruch 1984; CassatcUa 
& Gonzalez-Riestra 1990), while for V841 Oph it 
is E{B -V) = 0.30-0.58 (Bruch 1984; Cassatclla 
& Gonzalez-Riestra 1990; Weight et al. 1994; Vcr- 
bunt et al. 1997). Nova shells were not detected in 
Ha for either CV (Cohen 1985). In X-rays, DI Lac 
was undetected in the ROSAT PSPC All Sky Sur- 
vey, with a 2(7 upper limit of 0.013 cts in the 
11-201 channel (0.1-2.0 keV) energy range during 
299 s total exposure time, while V841 Oph was 
a marginal detection with 0.017 ± 0.007 cts 
in the "hard" 0.9-2.0 keV ROSAT energy range 



(0±0.025 cts s-i in the full ROSATcncigy range) 
during 402 s total exposure time (Verbunt et al. 
1997). 

We report here on the results of ground-based 
optical photometry and spectroscopy, and Rossi 
X-ray Timing Explorer [RXTE) X-ray observa- 
tions of DI Lac and V841 Oph. The optical pho- 
tometry comprises an almost decade-long light 
curve for each star, while the contemporaneous 
spectroscopy and X-ray observation repeatedly 
sampled each nova during separate intervals of 
« 45-55 d in length. 

2. Observations and Analysis 

Our long-term optical photometry of DI Lac 
and V841 Oph was accomphshed from 1990-1998, 
while the spectroscopic amd X-ray observations 
took place in 1997. These observations are sum- 
marized in Tables 1 and 2, and are discussed in 
detail below. 

2.1. Optical Photometry 

Our optical photometry data were acquired by 
RoboScope, a 41-cm telescope in Indiana equipped 
for automated differential CCD stellar photom- 
etry (Honeycutt & Turner 1992). All observa- 
tory operations (including data reductions) are ac- 
complished as fully unattended and unsupervised 
tasks, which makes practical the acquisition of 
long homogeneous data streams. Typically, Ro- 
boScope obtains one or two 4-min exposures per 
clear night for each of ^ 140 program stars. The 
data are reduced using the method of incomplete 
ensemble photometry (Honeycutt 1992). For DI 
Lac, 85 ensemble stars were used in 987 exposures 
over 9 observing seasons from 1990 November to 
1998 November. For V841 Oph, 24 ensemble stars 
were used in 611 exposures over 8 observing sea- 
sons from 1991 May to 1998 September. The zero- 
points of the differential light curves were estab- 
lished using secondary standard stars from Hen- 
den & Honeycutt (1995, 1997). Six such secondary 
standards were used for V841 Oph, while 13 were 
employed for DI Lac. Typical Icr uncertainties 
in the calibrated magnitudes are on the order of 
0.01-0.05 mag. Over the entire multiyear range of 
the RoboScope light curves, the mean magnitudes 
were V « 14.55 for DI Lac and V « 13.6 for V841 
Oph, with full ranges of variability of AF « 0.9 
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Fig. 1. — The optical (V) light curve of DI Lac 
from RoboScope. The top panel shows the en- 
tire light curve, while the bottom panel shows the 
region around the RXTE observations. The mid- 
point time of the RXTE observations is marked 
with an arrow. The horizontal bars above the 
data in the top panel show the eight observing 
seasons (the short initial observing season around 
HJD 2448200 is not marked), while the bar below 
the data shows the range of HJD displayed in the 
bottom panel. 

mag and « 1 mag, respectively. The Robo- 
Scope light curves of DI Lac and V841 Oph are 

shown in Figures 1 and 2. 

With only limited excursions to anomalously 
fainter magnitudes, the typical minimum magni- 
tude of « 14.7 for DI Lac has a sharp boimdary at 
any given epoch. However, the minimum bright- 
ness displays a slight, apparently linear, trend to- 
wards increasing brightness over the 8 years of Ro- 
boScope observation. If we ignore for the moment 
the small number of data points that fall below 
the well-defined low brightness limit seen in the 
top panel of Figure 1, then the typical minimum 
brightness of DI Lac changes from V « 14.75 mag 
at the beginning of the RoboScope coverage, to 
V fs 14.67 mag at the end. This corresponds 



Fig. 2.— As in Figure 1, but for V841 Oph. 

to a mean increase of the minimum brightness by 
AV » —0.01 mag yr~^. The maximum brightness 
of DI Lac is much less uniform, and it is difficult to 
determine if the maximum brightness behaves in 
the same manner as the minimum brightness. A 
linear fit to the complete DI Lac light curve (with 
each data point weighted by the inverse square of 
its Icr uncertainty) gives the rate of change of the 
mean magnitude as < -1-0.001 mag yr~^; that is, 
the mean brightness of DI Lac is essentially con- 
stant. However, if we exclude the 2 excursions 
to fainter magnitudes at HJD 2450600 and HJD 
2451000 (as well as the excursion to brighter mag- 
nitudes at HJD 2450300), then the slope of the 
linear fit changes to AV = —0.004 mag yr~^; that 
is, a trend towards increasing mean brightness. 

The faint magnitude excursions in the light 
curve of DI Lac are very similar to the "dips" seen 
by Honcycutt et al. (1998a) in long-term Robo- 
Scope light curves of five old novae and novalike 
CVs. In the latter systems, the dips were found to 
often be paired with a preceding or following out- 
burst, although the dips also sometimes occurred 
as isolated events. Honcycutt et al. (1998a) do 
not find a clear mechanism (e.g. disk instabilities, 
truncated disks, mass transfer modulation) that 
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is responsible for the dips, but conclude that the 
overall photometric behavior (dips and outbursts) 
of the novae and novalikes that they studied is 
likely to be governed by either a combination of 
disk and mass transfer events, or by mass transfer 
events alone. 

In contrast to the linear trend in the minimum 
brightness of DI Lac, the minimum brightness of 
V841 Oph appears to follow an almost sinusoidal 
trend that completes approximately 1.5-2 cycles 
during the 7.5 year RoboScope coverage. As seen 
in the top panel of Figure 2, the minimum bright- 
ness of V841 Oph varies from V « 14.5 mag to 
V w 13.8 mag. The maximimi brightness of V841 
Oph, which ranges from V » 13.1 mag toV ~ 13.4 
mag, displays the same behavior; that is, when the 
minimum brightness at a given epoch is at its faint 
(bright) level, the maximum brightness is also at 
its faint (bright) level. We note that the magni- 
tude range spanned by the minimum brightness in 
V841 Oph is about twice as large as that spanned 
by the maximum brightness. Although a linear fit 
to this overall trend is perhaps not the optimum 
choice, the corresponding rate of mean magnitude 
change is AV = +0.010 mag yr~^. It is possible 
(but not required by the data!) that the appar- 
ently linear trend in the light curve of DI Lac may 
also be a sinusoid, but with a much longer cycle 
length than in V841 Oph. 

We searched for periodicities in the range 10- 
100 days in the RoboScope light curve data in 
two ways: by applying the phase dispersion mini- 
mization (PDM) algorithm (Stellingwerf 1978), as 
well as performing an independent power spec- 
trum analysis using the CLEAN algorithm (Roberts 
et al. 1987). We performed the period search on 
both the entire data set for each CV, as well as 
on the data subsets from individual observing sea- 
sons (indicated in Figures 1 and 2). The most 
prominent period detected in each data set is listed 
in Table 3 (in several cases, most notably for the 
combined data sets for each nova, two significant 
periods are listed). In general, the two period 
search methods gave equivalent results to within 
better than 1-4 days. All of the detected time 
scales appear to correspond to quasiperiodic be- 
havior rather than truly periodic variability such 
as that resulting from, for example, an eclipsing 
orbit. The time scales have significant "flexibil- 
ity" (on the order of 1-5 days for periods < 50 d, 



5-20 days for periods > 50 d) over the multi-year 
length of the light curves (as shown by substan- 
tially broadened dips and peaks in the PDM and 
CLEAN analyses, respectively). Phase-binned light 
curves of both CVs folded on the shorter period 
found for each combined data set are shown in Fig- 
ure 3. We note that the true amplitude of varia- 
tion in a given cycle or observing season will be 
larger than the 0.1-0.2 mag suggested by Figure 
3, in which data from many cycles with slightly 
different amplitudes and periods have been av- 
eraged together. Inspection of Figures 1 and 2 
shows amplitudes of 0.4 0.8 mag for the full range 
of variability within each observing season. The 
long-term sinusoidal trend in the light curve of 
V841 Oph is also well-fit by a periodicity in the 
range 1800-1900 d or, with slightly less agree- 
ment, 1250-1300 d. (We note that these two pe- 
riod ranges are approximately related by the ratio 
3:2, so we are likely seeing an aliasing effect due to 
the fact that only 1.5-2 cycles of these long periods 
are contained in the light curves.) 

2.2. Optical Spectroscopy 

We obtained optical spectra of both CVs using 

the Double Imaging Spectrograph on the Apache 
Point Observatory (APO) 3.5-m telescope (Gille- 
spie et al. 1995)^ during University of Washington 
share time. Spectra of DI Lac were obtained in 
July and August of 1997; spectra of V841 Oph 
were obtained in May of 1997. The APO spectra 
have a resolution of « 2A and cover simultaneous 
wavelength ranges of 4200-5000A and « 5800- 
68OOA. The raw spectrum images were reduced in 
the standard fashion using IRAF. The instrumen- 
tal response was removed via spectra of standard 
stars (Massey et al. 1988) obtained on the same 
nights; however, slit losses due to guiding errors 
rendered the absolute flux calibration unreliable. 
Consequently, we have normalized the spectra to 
a constant continuum level of 1.0. 

A representative spectrum of DI Lac is shown in 
Figure 4. Note the prominent Balmer absorption 
troughs containing narrow, central emission cores. 
Our other two APO spectra of DI Lac (not shown) 
are essentially identical to this one, with the ex- 
ception that the Balmer emission cores in the Au- 
gust spectrum are somewhat stronger (relative to 

'also see http://www.apo.nmsu.edu/ 
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Fig. 3. — The RoboScopc light curves of DI Lac 
(top panel) and V841 Oph (bottom panel) folded 
on periods of 37 d and 36 d, respectively. A best- 
fit linear trend was subtracted prior to folding. 
The data have been averaged into phase bins of 
width 0.1, and are repeated twice. The horizontal 
error bars show the width of the phase bins, while 
the vertical error bars span ±3(Tn, where Cn is the 
standard deviation of the mean in each phase bin. 

He II) than in the July spectra. The He I A4471 
feature (see inset) is present with the same profile 
shape in all three of our spectra. It is suggestive 
of the absorption trough with central emission as 
seen in the Balmer lines. The feature labeled with 
"?" is likely an instrumental or reduction artifact 
since its wavelength does not match any identifi- 
able line and it occurs in all of our APO spectra 
(those of V841 Oph included). 

Additional optical spectra of both novae were 
obtained using the Hydra Multi-fiber Positioner -|- 
Bench Spectrograph (Harden et al. 1993, 1994) on 
the WIYN 3.5-m telescope^. Spectra of DI Lac 
were obtained in 1997 September, and spectra of 
V841 Oph were obtained in 1997 May. The WIYN 
spectra have a resolution of ~ lA and cover (us- 
able) wavelength ranges of w 6200-6800A for DI 
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Fig. 4. — Representative optical spectrum of DI 
Lac from APO on 1997 Jul 25 09:01:47 UT (HJD 
2450654.87744). The spectrum has been normal- 
ized to a continuum level of 1.0. The wavelengths 
of prominent line transitions are marked. The box 
inset in the left panel shows an expanded view of 
the region from 4390A to 4520A. The two features 
highlighted by dark underlines are marked (from 
left to right) "?" and "Hel" in the full spectrum 
plot. 

Lac and » 4400-5000A for V841 Oph. For the 
WIYN spectra of DI Lac, no standard star obser- 
vations were available, so it was not possible to 
correct for the instrumental response or perform 
a flux calibration. We do not show these spectra 
here, but note that they are qualitatively similar 
to the red APO spectrum shown in Figure 4 no 
differences corresponding to the different bright- 
ness levels in the RoboScope light curve are ap- 
parent. The V841 Oph WIYN data were reduced 
and calibrated as with the APO spectra. Again, 
slit losses and, additionally, the presence of cirrus 
made the absolute flux calibration unreliable, so 
we have normalized the spectra to a constant con- 
tinuum level of 1.0. All of the spectra of V841 Oph 
are shown in Figure 5 (the three WIYN spectra 
have been averaged together to improve the S/N). 
There is no indication of the absorption troughs 
seen in the spectra of DI Lac. Additionally, the 
C III/N III + He II emission complex appears to 
be stronger, with broader lines, in V841 Oph than 
in DI Lac. 
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Fig. 5. — Optical spectra of V841 Oph on (bottom 
to top): 1997 May 07 UT (APO), 1997 May 15 UT 
(APO), 1997 May 17 UT (WIYN), and 1997 May 
27 UT (APO). The spectra have been normahzed 
to a continuum level of 1.0. 



Flux-calibrated optical spectra of both DI Lac 
and V841 Oph are shown in Williams (1983), with 
resolutions of 5.lA and 6.7A, respectively, and in 
Ringwald et al. (1996), with resolution of 15-20A 
for both. Williams (1983) obtained his spectra 
in 1981, while Ringwald ct al. (1996) obtained 
theirs in 1990. In both papers, the spectra of 
both novae show steep blue continua; the observed 
continuum of DI Lac is bluer than that of V841 
Oph (but different reddening values could elimi- 
nate or reverse this characteristic). Ringwald et 
al. (1996) fit a power law of the form F\ oc A~" 
to the continuum regions 5600-6500A and 6650- 
8OOOA in their spectra. These fits yielded indices 
of a = 1.81(2) - dereddcned to = 3.10(3) using 
E{B -V) = 0.41 - for DI Lac and a = 1.45(1) - 
dereddened to ao = 2.79(1) using E{B-V) = 0.58 
- for V841 Oph. CassatcUa ct al. (1989) performed 
this power law fit to the wavelength region 1200- 
32OOA in lUE UV spectra of DI Lac and V841 
Oph. Using reddenings at the low ends of the 
ranges given in §1, E[B — V) = 0.15 for DI Lac 
and E{B -V) = 0.30 for V841 Oph, they calcu- 
lated dereddened slopes of ag = 1-5 and ao = 2.0, 
respectively. Comparison with the slope expected 
for the continuum flux distribution of a steady- 
state accretion disk composed of parcels radiating 



as blackbodics. a. = 2.3 (Pringlc 1981), suggests 
that the true reddening values for DI Lac and V841 
Oph are near the middle of the ranges given in §1. 

In 1936, both of these novae were observed spec- 
troscopically by Humason (1938), who noted their 
strong blue continua, but also noted the absence 
of any lines, emission or absorption. Since that 
time, the available published spectra of DI Lac and 
V841 Oph in Greenstein (1960), Williams (1983), 
and Ringwald et al. (1996) are consistent with our 
current spectra (see spectroscopic histories sum- 
marized in Ringwald et al. 1996), with the excep- 
tion that Williams (1983) did not report the de- 
tection of any He II A4686 emission in either star. 

2.3. X-ray Observations 

DI Lac and V841 Oph were each observed ten 
times during .luly-September and April-June of 
1997, respectively, with the Proportional Counter 
Array (PCA: Jahoda et al. 1996) on the RXTE 
sateUite (e.g. Bradt et al. 1993)^. The PCA con- 
sists of five xenon-methane proportional counters 
effective over the range 2-60 kcV (with 18% en- 
ergy resolution at 6 keV). Each RXTE visit con- 
stituted a Good-Time-Interval on-source lasting 
« 1000 2500 s. We extracted the net average X- 
ray count rate from each visit and the X-ray spec- 
trum from the ten combined visits to each star 
using the FTOOLS/LHEASOFT (v5.0) software 
package. We utilized the faint source background 
models (pca_bkgd_faintl7_e3vl9990824.mdl and 
pca_bkgd_faint240.e3vl9990909.mdl) and an up- 
to-date South Atlantic Anomaly history file 
(pca_saa_history_20000223). Extraction of the 
count rates and spectra was carried out as de- 
scribed in "The RXTE Cookbook" 9. 

The X-ray light curves constructed from the 
mean count rates at each visit for DI Lac and 
V841 Oph are shown in Figures 6 and 7, respec- 
tively, along with the corresponding sections of 
their RoboScope V band light curves. The light 
curve counts were summed in the energy range 
0-15 keV. Inspection of the X-ray spectra (see be- 
low) revealed an essentially zero count rate dom- 
inated by noise at energies higher than 15 keV. 
The mean X-ray count rates averaged over all ten 
visits to each target are 1.30 ± 0.32 ct s~^ for DI 



^also see http://heascirc.gsfc.nasa.gov/docs/xte/ 
^http://heasaj-c.gsfc.nasa.gov/docs/xte/recipes/cook_book.html 
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Fig. 6. — The optical {V) light curve of DI Lac 
from RoboScopc (small circles, left axis scale) with 
the mean RXTE count rate (0 15 kcV) at each of 
the ten epochs of observation (large circles con- 
nected by a dotted line, right axis scale). One- 
sigma error bars are shown on all of the points; the 
arrows mark the times when optical spectra were 
obtained. The lowered arrows mark the WIYN 
spectra; the others are from APO. 

Lac and 0.56 ± 0.73 ct s"! for V841 Oph. Con- 
sidering these very low count rates, it is not clear 
whether or not the variability seen in the X-ray 
light curves is real. 

Additional analysis and model-fitting (DI Lac 
only) of the X-ray spectra of these novae was per- 
formed using the routine XSPEC (vll). Although 
we extracted the full range of available energy 
channels in the X-ray spectra, only the range 0-15 
keV was used for fitting the models to the DI Lac 
data. The X-ray spectra for DI Lac and V841 Oph 
are shown in Figures 8 and 9. The X-ray spec- 
trum of V841 Oph is shown on the same scale as 
that of DI Lac; the former is clearly much weaker 
than the latter. We attempted to extract a bet- 
ter spectrum of V841 Oph by using only the data 
from the four visits with the largest count rates 
(visits 3, 5, 9, and 11 - see Table 2). Unfortu- 
nately, this spectrum was virtually identical to the 
total combined spectrum (i.e. flat and featureless), 



Fig. 7. — As in Figure 6, but for V841 Oph. 

with the exception of a small upward offset of the 
mean count rate, from about 0.02 cts s~^ keV~^ in 
the total combined spectrum to about 0.04 cts s^^ 
keV~^ in the "high" count rate spectrum. Conse- 
quently, no attempt was made to fit a model to 
the spectrum of V841 Oph. 

For DI Lac, we fit three simple models to the X- 
ray spectrum: (1) blackbody, (2) Raymond-Smith 
thermal plasma, and (3) brcmsstrahhmg. Each 
model was modified by a multiplicative component 
representing photoelectric (H column) absorption. 
More complex models were not warranted owing to 
the small number of counts per energy channel in 
the data. In all cases, the data were weighted us- 
ing the XSPEC-recommended scheme appropriate 
for small count numbers, Wi = l + {Ni-\- 0.75)°-^ 
(Gehrels 1986). The results of these model fits are 
summarized in the first three rows of Table 4. All 
of the models produce nearly indistinguishable fits 
to the observed spectrum. The bremsstrahlung 
model is formally the best-fitting (with lowest re- 
duced and highest null hypothesis probability) 
and is plotted over the observed spectrum in Fig- 
ure 8. 

All of the models were only able to constrain 
the hydrogen column density to an upper limit; 
the values of TVh quoted in the first three rows of 
Table 4 are the upper ends of the 1-parameter 90% 
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Fig. 9. — The combined X-ray spectrum of V841 
Oph on the same scale as the X-ray spectrum of 
DI Lac shown in Figure 8. 



Fig. 8. — Top panel shows the combined X-ray 
spectrum of DI Lac (crosses) with the best-fitting 
bremsstrahlung model modified by photoelectric 
(H column) absorption (solid line). Bottom panel 
shows the residuals to the model fit. 

confidence intervals. We used the relation 

Nn = E{B -V)x 5.8 x lO^^cm-^mag"^ 

(Bohlin et al. 1978) to estimate values of A^h span- 
ning the reddening range given in §1. Wc then 
re-fit the models to the X-ray spectrum of DI Lac 
with only two free parameters: temperature and 
normalization. The results arc listed in the bot- 
tom six rows of Table 4. While slightly different 
from those obtained from the original model fits, 
they arc completely consistent with them. The 
bremsstrahlung model yields a slightly lower nomi- 
nal temperature of fcT w 4.1-4.2 keV (vs. kT = 4.5 
keV when A^h is a free parameter) and is still the 
best- fitting model for both values of A^h- 

We note that the uncertainties quoted in Table 
4 for kT and Cnorm arc 1-paramctcr 90% confi- 
dence intervals; the interdependent parameter un- 
certainties will be slightly different. For example, 
for the bremsstrahlung model (with free A^h), the 
XSPEC routine steppar^° gives a 2-parameter 90% 



confidence interval of 

kT = 4.5t J:^ keV 
for kT with respect to A^h, and 

kT = 4.5tJ;^ keV 

for kT with respect to Cnorm- 

The X-ray cormt rates were too low to reliably 
fit spectrum models to any of the ten individual 
RXTE visits to DI Lac. We constructed "low" 
and "high" state X-ray spectra by combining the 
data from visits with the lowest mean count rates 
(1, 2, 3, 4, 7) and from visits with the highest 
mean count rates (5, 6, 9) to see if there were any 
differences between the low and high X-ray bright- 
ness states^^. Other than an « 1 cts s~^ keV~^ 
increase at the peak energy of the high state spec- 
trum, the two spectra are essentially indistinguish- 
able from each other (and from the combined spec- 
trum from all ten visits). We conclude from this 
that either the spectral energy distribution of the 
X-ray emitting source in DI Lac does not change as 
the overall brightness of the X-ray source changes, 



^"This routine iteratively refits the model to the data while 



stepping the values of selected parameters through a given 
range and produces a plot in parameter space of the 2- 
dimensional contours. 
^^Visits 8 and 10 were not used. The former was deemed 
to be an "intermediate ' brightness state, and the latter an 
anomalously low state. 
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or that the variability seen in the X-ray hght curve 
of DI Lac is not real (i.e. it is a side-effect of the 
very low count rates). 

3. Discussion 

3.1. Optical Photometric Variability 

Inspection of the RoboScopc light curves of DI 
Lac and V841 Oph reveals variability on three dis- 
tinct time scales. First, in the range 10 100 d, 
quasiperiods with typical lengths of ^ 30 50 d are 
present in each of these CVs. Yet, the presence (or 
lack) of periodic and/or quasiperiodic variability 
on time scales up to a few hundred days in V841 
Oph has been contested in the literature. (To 
the best of our knowledge, no extensive database 
of photometric observations of DI Lac has been 
published prior to this work.) Using 420 archival 
visual observations of V841 Oph spanning 28 yr, 
Delia Valle & Rosino (1987) quote an average pe- 
riod of 51.5 d with individual cycles ranging from 
45 d to 57 d. Shara et al. (1989) observed V841 
Oph once per night for 31 nights during a 45 night 
interval in 1986 using a CCD camera + B filter, 
and found a period of ^ 40 d. All of these values 
are comparable to the ^ 35-50 d we found for our 
complete RoboScope light curve of V841 Oph, as 
well as to the range of periods found for individual 
observing seasons in the RoboScope data (see Ta- 
ble 3). On the other hand, Duerbeck (1992) notes 
the presence of considerable variability in archival 
light curves of V841 Oph (including some of the 
same data used by Delia Valle & Rosino 1987), but 
finds no evidence for periodicities up to 200 d. 

It is quite interesting that the variability time 
scales in these two novae are similar to those that 
characterize dwarf nova outbursts (Warner 1995, 
ch. 3). Old novae are presumed to have mass 
transfer rates, M, above the threshold level (Merit) 
at which the accretion disk thermal instability 
mechanism operates in dwarf novae to produce 
their outbursts (Osaki 1996). However, the hiber- 
nation theory of cyclic evolution between CV vari- 
ability types (Shara 1989, and references therein), 
predicts that when M eventually decreases back 
into the instability regime some decades after the 
nova outburst, then dwarf nova outbursts should 
resume. Both GK Persei (= Nova Per 1901; Bian- 
chini et al. 1986) and V446 Herculis (= Nova Her 
1960; Honeycutt et al. 1998b) have displayed dwarf 



nova eruptions in their post-nova stages, and sim- 
ilar behavior has been suggested for a number 
of other old novae, commencing « 50-200 yr af- 
ter their outbursts (Livio 1989; Warner 1995, ch. 
4). The amplitudes of variability displayed by 
DI Lac and V841 Oph (see Figures 1 and 2) are 
smaller than those typically observed for dwarf no- 
vae {AV ~ 2-4 mag), but we cannot rule out the 
possibility of a disk instability mechanism being 
in operation after all, it has been 90 years since 
the nova outburst of DI Lac and over 150 years 
since that of V841 Oph. Further, the classic ther- 
mal instability mechanism is not well-explored for 
high-M disks, and other disk instabilities might 
operate as well (e.g. Whitehurst 1988; Mineshige 
1993; Godon 1998). 

The amplitudes of the ~ 30-50 d oscillations 
in DI Lac and V841 Oph are similar to those of 
the spaced stunted outbursts seen in some nova- 
like CVs (Honeycutt et al. 1998a). One of the sug- 
gestions for stunted outbursts made in Honeycutt 
et al. (1998a) is that they are dwarf nova erup- 
tions seen against a background of brighter light 
in the system. If that suggestion is correct, then 
the oscillations reported here (as well as similar 25 
d oscillations reported for RW Tri, Honeycutt et 
al. 1994) may be examples of dwarf nova type out- 
burst behavior at relatively short recurrence times, 
Tn < 50 d (Tn is typically > 10 d to hundreds of 
days - even tens of years in extreme cases - for 
dwarf novae; Warner 1995, ch. 3). Among the true 
dwarf novae, the Z Camelopardalis stars have out- 
burst recurrence times of « 10-30 d, somewhat 
shorter than the time scales for the oscillations 
observed in these novae. The Z Cam stars are 
thought to have high mass transfer rates (near 
Merit)- This suggests that if the oscillations ob- 
served in novae are analogous to dwarf nova out- 
bursts, then they may be linked to the presence 
of relatively high mass transfer rates that are, 
nonetheless, smaller than both Merit and M in the 
Z Cam stars. 

Second, V841 Oph displays a distinct sinusoidal 
variation with a best period of 1800-1900 d (4.9- 
5.2 yr) or a slightly less preferred period of 1250- 
1300 d (3.4-3.6 yr). DI Lac does not appear 
to display similar sinusoidal variability, unless it 
occurs with a much longer period than in V841 
Oph. Bianchini (1990) reports a long-term period 
in V841 Oph of 3.4 yr with mean amplitude of 



9 



^ 0.3 mag, which was determined from 420 visual 
observations spanning 30 yr (essentially the same 
early- to mid-20th century data utilized by Duer- 
beck 1992, see below). Richman ct al. (1994) re- 
analyzed the light curves presented by Bianchini 
(1990) and consider the case for a multi-year pe- 
riodicity in V841 Oph to be very weak. However, 
with our RoboScope light curve, sinusoidal vari- 
ability on comparable multi-year times scales has 
now been observed in two separate data sets for 
V841 Oph spanning almost 80 years of observa- 
tion. While this still does not provide firm ev- 
idence of strictly periodic behavior, it certainly 
points to the presence of a mechanism operating in 
this CV that modulates its brightness on a char- 
acteristic time scale of several years. 

Photometric variability in old novae on time 
scales of several years has generally been at- 
tributed to solar-type magnetic cycles on the sec- 
ondary star that might be able to control the rate 
of mass loss through the inner Lagrangian point 
(Bianchini 1990; Applegate 1992). Although Rich- 
man et al. (1994) find that no CV studied to date 
displays strictly periodic behavior over multi-year 
intervals, they also note that the observed ampli- 
tudes (~ 0.2 mag) and "apparent" time scales of 
variability of 5-40 yr are "plausible consequences 
from solar- type magnetic cycles." We note that 
the shorter time scale (~30-50 d) variability dis- 
cussed above might also be linked to modulation 
of M due to secondary star magnetic activity. If 
solar-type magnetic cycles can affect M on time 
scales of years, then starspots induced by the mag- 
netic cycle that migrate under the Li point might 
also affect M on shorter time scales, in a manner 
suggested by Livio & Pringle (1994) to explain the 
very low brightness states seen in some CVs. 

Third, the mean magnitudes of both CVs are 
changing slowly, at rates measured in a few milli- 
magnitudes (mmag) per year. The mean bright- 
ness of DI Lac is increasing by 4 mmag yr^^ since 
1990. Duerbeck (1992) analyzed a large set of vi- 
sual, photoelectric, and CCD observations of post- 
novae available in the literature (see references in 
Duerbeck 1992). These data cover a large part 
of the 20th century, from the 1920's through the 
1980's. Duerbeck (1992) summarizes the long- 
term behavior of DI Lac as exhibiting a steady 
decline of 13 ± 1 mmag yr^^ in visual observations 
obtained between 1921 and 1952. Sparse photo- 



electric and CCD data from 1953 to 1981 sug- 
gest a possible small brightness increase, which 
is supported by our recent RoboScope observa- 
tions. For V841 Oph, visual observations from 
1919-52 and 1978-91 yield a brightness decline of 
7± 1 mmag yr~^, while "no definitive brightness 
decline" (Duerbeck 1992) was found in photoelec- 
tric data from 1954 to 1988. Our post-1991 ob- 
servations of V841 Oph suggest a trend towards 
decreasing mean brightness comparable to that 
noted by Duerbeck (1992) in the archival visual 
observations, although the obvious sinusoidal vari- 
ation that dominates our light curve at long time 
scales makes this conclusion somewhat suspect. 

In the "classical" hibernation scenario (Shara 
1989), novae stay bright for ~ 50-300 yr following 
outburst, gradually declining in brightness dur- 
ing this time to a hibernating state character- 
ized by low brightness (My ~ 10 or fainter) and 
low mass accretion rate {M < 10~^^ Mq yr~^). 
However, the photometric record of DI Lac from 
the 1920's (a decade after its outburst) to the 
1990's does not show a simple decline in bright- 
ness. Instead, the brightness of DI Lac decreased 
for about four decades after outburst, but then ap- 
parently leveled off for several decades. In the last 
decade (covered in our RoboScope light curve), 
the brightness of DI Lac has been increasing. It is 
possible that we are seeing the influence of a mech- 
anism (such as magnetic activity of the secondary 
star - see above - and/or accretion-induced ir- 
radiation of the secondary star) that modulates 
M (and, hence, the brightness) in DI Lac over 
a long-term cycle (with period on the order of 
many decades) that obscures the general decline 
predicted by the hibernation scenario. 

3.2. Optical Spectrum Variability 

The optical spectra of DI Lac and V841 Oph 
are quite different, and this can possibly be as- 
cribed to their difference in post-outburst ages. 
The spectrum of V841 Oph is similar to those of 
novalike CVs (e.g. Warner 1995, ch. 4) - this sug- 
gests that V841 Oph still has a high accretion rate, 
even « 150 yr after outburst. The narrow emis- 
sion components in DI Lac are suggestive of emis- 
sion originating from the irradiated inner face of 
the secondary star, while the absorption troughs 
imply the presence of optically thick material in 
the system. The latter feature could be material 
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ejected during the outburst of this younger post- 
nova; the lack of any detected Ha emission shell 
(Cohen 1985) does not preclude the existence of 
denser, non-cniitting circumstcUar material. The 
former feature offers a possible explanation for the 
gradual brightness increase seen in our RoboScope 
light curve of DI Lac (discussed in §3.1) if the sec- 
ondary star is being slowly heated via irradiation 
and is, in turn, increasing the rate of mass transfer 
through the Li point. 

Although the RoboScope coverage is incom- 
plete, we infer from the adjacent data that DI 
Lac was faint during our July spectra and bright 
during our August spectrum (see Figure 6). As 
mentioned in §2.2, the only difference between our 
July and August spectra is that the emission cores 
were somewhat stronger in August (bright state) 
than July (faint state). This is consistent with 
the hypothesis that these narrow emission features 
originate on the irradiated face of the secondary 
star if we make the logical assumption that the 
irradiation increases when DI Lac is bright. It 
does not, however, illuminate the exact mecha- 
nism producing the irradiation (i.e. whether the 
"excess" flux in the high state originates near the 
disk center - presumably due to a disk instability 
producing increased accretion onto the WD - or 
in the outer disk - presumably due to an increase 
in mass transfer through the Li point.) If we as- 
sume that the spectrum shown in Figure 4 (from 
July 27 09:01:47 UT) was obtained at an arbitrary 
orbital phase of 0.0, then the August 17 spectrum 
was obtained at a relative orbital phase of 0.9 (us- 
ing Porb = 0.543773 d; Ritter & Kolb 1998). Thus, 
these two spectra were obtained at similar orbital 
phases, and we do not expect the difference in the 
narrow emission component strength to be only 
due to system orientation. 

The RoboScope coverage of V841 Oph during 
our spectroscopic (and X-ray) observations is more 
sparse than for DI Lac, but we can infer that the 
CV was returning to the faint state during our 
May 07 spectrum, was in its faint state during 
our May 14 and 17 spectra, and was likely near 
the bright state during our May 27 spectrum (sec 
Figure 7). (The two elevated brightness points at 
HJD 2450610 and HJD 2450625 suggest that V841 
Oph may have returned to a bright state sometime 
between the low states bracketing HJD 2450590 
and HJD 2450630.) If this is the case, then the 



somewhat stronger Ha emission of the May 07 and 
May 27 spectra (see Figure 5) may be linked to 
the bright state of the CV. As with DI Lac, we 
calculated relative orbital phases for each of these 
spectra (using Porb = 0.60423 d; Ritter & Kolb 
1998), and obtained (j) = 0.0, 0.6, 0.75, 0.1 for May 
07, 14, 17, 27, respectively. Unfortimately, this 
casts some doubt on the link between Ha emis- 
sion strength and brightness state, since both of 
the bright state spectra (phases 0.0, 0.1) were ob- 
tained at different orbital phases than the faint 
state spectra (phases 0.6, 0.75). So, we cannot 
rule out the influence of system orientation effects 
in these spectra. 

3.3. X-ray Variability 

Unfortunately, because of the weak X-ray emis- 
sion from these CVs, little can be firmly stated 
about their X-ray variability. The count rates dur- 
ing visits 5 and 6 to DI Lac, which took place when 
we infer from the RoboScope light curve that the 
CV was in a bright state, are slightly elevated com- 
pared to the preceding visits. However, visit 7 
also occurred during this optical bright state and 
does not show an elevated count rate. The mean 
X-ray count rate during each visit to V841 Oph 
also does not display any strong correlation with 
the corresponding optical state. The lack of large 
changes in X-ray flux or spectrum during the op- 
tical variations (as usually evident in dwarf novae; 
e.g. Szkody et al. 1999), argues against a disk in- 
stability scenario. Alternatively, this could indi- 
cate that the optical brightness is determined by 
activity in the outer disk only and, therefore, is 
not reflected in the X-ray behavior; however, the 
X-ray count rates are too low (and their corre- 
sponding error bars too large) to make any firm 
conclusions. 

4. Conclusions 

Our long-term optical light curves of the novae 

DI Lac and V841 Oph obtained with RoboScope 
reveal quasiperiodic variability with a characteris- 
tic time scale of ~ 30-50 d in both CVs. In addi- 
tion, the light curve of V841 Oph displays evidence 
for sinusoidal variability with a period of 3.5-5 yr. 
The latter cannot be said to be strictly periodic 
since our data set covers only « 1.5 cycles; how- 
ever, when this detection is added to past reports 
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of multi-year periodicities in V841 Oph that have 
been reported in the literature, a strong case can 
be made for the presence of repeating multi-year 
variability with a preferred time scale in this CV. 
The most likely origin of such behavior is a solar- 
type magnetic cycle of the secondary star that 
modulates the rate of mass transfer through the 
Li point. If this is the case, then the shorter ~ SO- 
SO d quasiperiodic variability in this system might 
also be related to the magnetic activity on the sec- 
ondary star; for example, due to starspots that mi- 
grate under the Li point and temporarily throttle 
mass transfer. DI Lac does not show evidence for 
cyclic multi-year variability, unless it occurs with 
a much longer period than the length of our Robo- 
Scope coverage (> 10 yr). This casts some doubt 
on the origin of the ~ 30-50 d quasiperiodic vari- 
ability as a facet of the secondary star's magnetic 
activity, since both V841 Oph and DI Lac show 
the ~ 30-50 d quasiperiodic behavior, but only 
V841 Oph displays evidence for a multi-year solar- 
type magnetic cycle on the secondary star. On 
the other hand, the ~ 30-50 d time scale is very 
reminiscent of that expected for dwarf nova type 
behavior (although the amplitude of variability in 
these novae is smaller than in typical dwarf nova 
outbursts). This raises the possibility that this 
variability is caused by a disk instability (either 
the thermal disk instability that leads to dwarf 
nova outbursts operating at a low level or some 
other form of disk instability). 

The X-ray spectrum of DI Lac is fit almost 
equally well by the three models we tried: a sim- 
ple blackbody, a Raymond-Smith thermal plasma, 
and bremsstrahlung emission. More complicated 
models (e.g. involving multiple components, lines, 
etc.) arc unwarranted due to the low X-ray flux. 
In addition to being the most physically plausi- 
ble X-ray emission emission mechanism in a (non- 
magnetic) CV, the bremsstrahhmg model is for- 
mally the best-fitting. Our X-ray spectrum model 
fit parameters (kT ~ 4 keV) are consistent with 
those obtained from ROSAT X-ray spectra of a 
sample of 37 disk-accreting CVs (Richman 1996). 
Unfortunately, the X-ray count rates in both of 
these systems were too low for any conclusions 
to be made about their time-resolved X-ray be- 
havior. X-ray observations sample the innermost 
disk region, and we would expect to see differ- 
ent time-resolved behavior as the novae go into 



their optically bright state if the transition is trig- 
gered by a disk instability vs. a change in M from 
the secondary star. Additional time-resolved X- 
ray observations using the more sensitive Chandra 
and/or XMM X-ray satellites may be necessary to 
illuminate the inner workings of these old novae. 
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Table 1 

Log of Observations for DI Lac 



UT Date and Time'' 


HJD'^ 


Observatory 


Mode^ 


Exposure 

or GTI (h) 


Comment 


1990 Nov 12 


2448207.7 


RoboScope 


P 


240^^ 


starting date 


1997 Jul 25 10:43:44 


2450654.94822 


RXTE 


X 


2528 


ID 20037-02-01-00 


1997 Jul 27 09:01:47 


2450654.87744 


APO 


s 


900 




1997 Jul 27 09:16:53 


2450654.88792 


APO 


s 


737 




1997 Jul 30 04:16:00 


2450659.67921 


RXTE 


X 


944 


ID 20037-02-02-00 


1997 Aug 03 04:11:12 


2450663.68424 


RXTE 


X 


1536 


ID 20037-02-03-00 


1997 Aug 07 04:12:00 


2450667.67683 


RXTE 


X 


1680 


ID 20037-02-04-00 


1997 Aug 12 03:42:40 


2450672.65668 


RXTE 


X 


592 


ID 20037-02-05-00 


1997 Aug 17 03:42:18 


2450677.65666 


APO 


s 


900 




1997 Aug 17 03:45:36 


2450677.65893 


RXTE 


X 


1024 


ID 20037-02-06-00 


1997 Aug 23 03:46:40 


2450683.65991 


RXTE 


X 


832 


ID 20037-02-07-00 


1997 Aug 27 03:46:56 


2450687.66024 


RXTE 


X 


736 


ID 20037-02-08-00 


1997 Sep 01 08:12:48 


2450692.84503 


RXTE 


X 


1616 


ID 20037-02-09-00 


1997 Sep 03 05:56:28 


2450694.74974 


WIYN 


s 


900 


No flux calibration 


1997 Sep 04 04:33:17 


2450695.69269 


WIYN 


s 


900 


No flux calibration 


1997 Sep 06 08:08:48 


2450697.84239 


RXTE 


X 


1776 


ID 20037-02-10-00 


1998 Nov 30 


2451147.6 


RoboScope 


p 


240^^ 


ending date 



At mid-point of exposure or Good-Timc-Intcrval (GTI) 
b p = ground-based photometry; S = ground-based spectroscopy; X = X-ray 
1-2 measurements night """^ 
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Table 2 

Log of Observations for V841 Oph 



UT Date and Time'* 


HJD'' 


Observatory 


Mode'' 


Exposure 
or GTI (s) 


Comment 


1991 May 31 


2448407.7 


RoboScope 


P 


240^= 


starting date 


1997 Apr 30 16:10:57 


2450569.17897 


RXTE 


X 


1936 


ID 20037-01-01-00 


1997 May 05 12:50:25 


2450574.03997 


RXTE 


X 


1856 


ID 20037-01-02-00 


1997 Mav 07 Ofi-30-32 


2450575 771 20 


APO 





600 




1997 May 10 20:59:13 


2450579.37966 


RXTE 


X 


1936 


ID 20037-01-03-00 


1997 May 14 06:37:42 


2450582.77618 


APO 


s 


600 




1997 May 15 21:12:33 


2450584.38911 


RXTE 


X 


1408 


ID 20037-01-04-00 


1997 May 17 09:35 


2450585.899 


WIYN 


s 


4500 


3 combined spectra 


1997 May 20 07:58:41 


2450588.83795 


RXTE 


X 


1792 


ID 20037-01-05-00 


1997 May 26 06:28:11 


2450594.77523 


RXTE 


X 


1872 


ID 20037-01-06-00 


1997 May 27 06:22:21 


2450595.76552 


APO 


s 


600 




1997 May 30 05:06:57 


2450598.71888 


RXTE 


X 


1648 


ID 20037-01-07-00 


1997 Jun 04 22:57:37 


2450604.46244 


RXTE 


X 


832 


ID 20037-01-08-00 


1997 Jun 09 06:43:45 


2450608.78613 


RXTE 


X 


2160 


ID 20037-01-09-00 


1997 Jun 18 16:24:33 


2450618.18935 


RXTE 


X 


1376 


ID 20037-01-11-00'^ 


1998 Sep 04 


2451060.5 


RoboScope 


p 


240^= 


ending date 



At mid-point of exposure or Good- Time-Interval (GTI) 
b p = ground-based photometry; S = ground-based spectroscopy; X = X-ray 
1-2 measurements night ~^ 

Our tenth RXTE observation of V841 Oph on 1997 Jun 14 UT (ID 20037-01-10-00) was unusable 
due to a software event that shut down all of the detectors during the slew to the target. 



Table 3 

Periods in Optical Light Curves 



Observing 
Season 


DI Lac 
Period (d) 


V841 Oph 
Period (d) 


1 


40 


35 


2 


29 


36 


3 


39 


69 


4 


31 


40 


5 


38 


42 


6 


75,40 


49 


7 


80,28 


54 


8 


31,21 


48 


combined 


43,37 


49,36 
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Table 4 

X-RAY Spectrum Models for DI Lac 



Model 


Reducedx^ 


Po^ H Column 2 15 kcV 


Parameters'^ 




(47 dof)^ 


Density, iVn Model Flux 


kT Cnorm 






(10^^ cm~^) (erg cm^^) 


(keV) 



Blackbody 





.705 


0.937 


< 1.2 


2, 


,2 


X 


10- 


-12 


1. 


1-1-0.1 
^-0.2 


3, 


0+0.6 
■^-0.2 


X 


10- 


-5 


Raymond-Smith'^ 





.617 


0.982 


< 3.3 


2. 


.6 


X 


10- 


-12 


3.: 


5-1-0.7 

^-1.1 


3. 


4+3.5 
■^-0.6 


X 


10- 


-3 


Bremsstrahlung 





.484 


0.999 


< 1.8 


2. 


.7 


X 


10- 


-12 


4.1 


c: -1-1.3 
3-1.3 


1, 


1+0.7 
•■■--1.1 


X 


10- 


-3 


H Column Density 


fixed 


from 


E{B - V) 


= 0.15: 


























Blackbody 





.683 


0.954 


0.09 


2, 


,2 


X 


10- 


-12 


1.: 


1-1-0.1 
J^-0.1 


3, 


1+0.4 
•-'^-0.3 


X 


10- 


-5 


Raymond-Smith'* 





.604 


0.986 


0.09 


2, 


,7 


X 


10- 


-12 


2.1 


q -1-0.9 


3, 


7+0.9 
•'-0.9 


X 


10- 


-3 


Bremsstrahlung 





.471 


0.999 


0.09 


2, 


.7 


X 


10- 


-12 


4.: 


5-1-1.5 
^-1.0 


1, 


9+0.4 
•^-0.3 


X 


10- 


-3 


H Column Density 


fixed 


from 


E{B - V) 


= 0.41: 


























Blackbody 





.689 


0.951 


0.24 


2, 


.3 


X 


10- 


-12 


1. 


1-1-0.1 
J--0.1 


3, 


q+0.3 
■"J-0.4 


X 


10- 


-5 


Raymond-Smith'* 





.606 


0.986 


0.24 


2, 


.7 


X 


10- 


-12 


2.! 


q-f-0.9 


3, 


O+0.9 
■^-0.9 


X 




-3 


Bremsstrahlung 





.475 


0.999 


0.24 


2, 


.7 


X 


10- 


-12 


4.: 


1-1-1.4 
J--1.0 


1, 


9+0.4 
•^-0.3 


X 


10- 


-3 



^ 48 dof for fixed A'h models. 

^ Null Hypothesis Probability: the probability of obtaining a value of greater than or equal to 
the observed if the model is correct. A value close to 1.0 indicates a good fit. 

Quoted uncertainties are 1-parametcr 90% confidence intervals. T is the plasma temperature 
for the bremsstrahlung and Raymond-Smith plasma models. In all models, Cnorm is a normalization 
constant. 

** Additional model parameters of redshift and abundances were fixed at 0.0 and the values from 
Feldman (1992), respectively. 
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